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MODIFICATIONS TO FOURIER’S LAW OF HEAT CONDUCTION 
THROUGH THE USE OF A SIMPLE COLLISION MODEL 

By Willard E. Meador 
Langley Research Center 

SUMMARY 

The exact first-order solution of an approximate kinetic equation is obtained for 
the problem of viscous gases flowing through circular pipes with constant temperature 
gradients. When expressed in terms of transport coefficients, the resulting flux-force 
relations are assumed to be correct since they agree in the proper limit with the Grad 
13 -moment approximation. However, the complete flux -force relations have important 
contributions from at least the first 22 moments; in particular, terms of convection mag- 
nitude are added to both the axial and radial components of the 13-moment equations for 
thermal conduction. These terms, which are linearly proportional to the magnitude of 
the flow velocity (instead of the square), significantly alter the conventional Navier- 
Stokes description of heat transfer in the problem considered. They also affect the 
evaluation of the pressure tensor. 


INTRODUCTION 

Calculations of heat transfer in flowing gases are normally based upon the conven- 
tional theory of Navier and Stokes (ref. 1), in which Newton’s law of viscosity and 
Fourier's law of thermal conduction are employed in the macroscopic equations of 
change. Despite the fact that the validity of such laws has been demonstrated in many 
practical applications and also on the basis of the second approximation to the kinetic 
theory of Chapman and Enskog (ref. 2), their general applicability to all near -equilibrium 
systems composed of ordinary monatomic gases is not guaranteed. The purpose of the 
present research to to investigate the accuracy with which the conventional theory 
describes a particular example of viscous pipe flow. 

Two generalizations of the simple Newton and Fourier laws are well known - 
namely, Grad's 13-moment addition of coupling terms (ref. 1, p. 495) to the heat flux 
and traceless pressure tensor. However, since the coupling terms can be obtained from 
the Chapman -Enskog third approximation, the assertion is frequently made that they are 
of little importance unless the gas properties vary appreciably in the distance of a mean 



free path. This argument implicitly assumes the mean free path to be a valid perturba- 
tion parameter and the divergence of the traceless pressure tensor to be negligible com- 
pared with the pressure gradient, even though simple examples can be cited to the 
contrary. 

One such example is the steady-state, pressure -driven, viscous flow of a gas 
through a pipe. The pressure gradient is balanced in this problem by the divergence of 
the traceless pressure tensor, which therefore cannot be regarded as a third- 
approximation force and should not be neglected in the velocity distribution function. In 
addition, the 13-moment coupling contribution to the conductive heat flux is directly pro- 
portional to the divergence of the traceless pressure tensor and is just as important as 
Fourier's term. Although neither contribution is significant when compared with heat 
transfer by convection, the question naturally arises as to whether higher Chapman- 
Enskog or higher Grad moment approximations yield conductive terms of much larger 
magnitude. The answer to this question is the primary objective of the present research. 

Since higher approximations to either the Chapman -Enskog or the Grad solutions 
of the Boltzmann integro-differential equation are very difficult to obtain, a common 
practice (see ref. 3, for example) is the substitution of simple collision models for the 
complete Boltzmann collision integrals in order to reduce the kinetic equations to differ- 
ential form. Perhaps the most frequently used model is that of Krook (ref. 3), which 
assumes the collision integrals to be directly proportional to the difference between the 
equilibrium and nonequilibrium velocity distribution functions. This model has several 
attractive features, but it does not distinguish between the effective collision frequencies 
for different transport phenomena; hence, the Krook predictions of quantities like the 
Prandtl number contain substantial errors. The model developed in the present research 
for obtaining a differential form of the kinetic equation avoids such difficulties. 

The first-order (that is, second-approximation) form of the resulting kinetic equa- 
tion is the fundamental relation used in the present study. It is obtained from the approx- 
imate model equation for viscous pipe flow by deleting all terms corresponding to squares 
and higher powers of the mean gas velocity and its driving force. This procedure is less 
restrictive than Chapman and Enskog's second approximation because spatial derivatives 
of the first-order perturbation function are not arbitrarily neglected in the first-order 
kinetic equation. Accordingly, the aforementioned 13 -moment additions to the traceless 
pressure tensor and the heat flux, as well as other phenomena of equal or greater impor- 
tance, are included automatically without reference to higher approximations or orders. 

The utility of the approximate first-order kinetic equation rests in the fact that 
exact solutions can be found for many problems - including the one of viscous gases 
flowing through circular pipes with constant temperature gradients. With the reasonable 
assumption that errors inherent in the collision model affect the values of effective 
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collision frequencies but not their ratios, the velocity moments of the exact solution are 
expected to give accurate flux-force relations if precise transport coefficients are used 
in the final expressions. Some support for this speculation is provided by the precise 
reduction to Grad's formulas in the 13 -moment limit and by the correct prediction of 
the Prandtl number. However, the complete flux -force relations have important contri- 
butions from at least the first 22 moments; in particular, contributions of the same mag- 
nitude as heat convection are added to both the axial and radial components of the 
13-moment equations for thermal conduction. These terms, which are proportional to 
the magnitude of the flow velocity (instead of the magnitude squared), significantly alter 
the conventional Navier -Stokes description of heat transfer in the problem considered. 
They also affect the evaluation of the pressure tensor. 

Previous efforts to obtain exact solutions of linearized Boltzmann equations have 
been made by Wang Chang and Uhlenbeck for the problems of heat transport between 
parallel plates (ref. 4), Couette flow (ref. 5), and flow near a surface (ref. 6). However, 
these solutions are not given in the closed form of the present research but appear 
instead as formal expansions in the complete set of spherical harmonics. Only the first 
few coefficients in these expansions are derived. Truesdell (ref. 7, see pp. 74-83) 
exactly solved a problem of simple shear in an infinite medium by considering the com- 
plete set of macroscopic equations of change and making assumptions about the behavior 
of the higher moments. These assumptions do not apply to the problem treated herein. 

SYMBOLS 


a radius of circular pipe 

ratio of collision integrals defined by equation (21) 
b impact parameter 

'c particle velocity relative to pipe 

Cp specific heat per particle at constant pressure 

e ij ,6 ij combinations of velocity or heat-flux derivatives (see eqs. (87) and (88)) 
f velocity distribution function 

f^ Maxwellian distribution function relative to mean gas flow 
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trial velocity functions 


s;s v ■ • -»g 9 

i,j,k indices; also used to denote vector and tensor components 

A A A 

i,j,k unit vectors along x-, y-, and z-axes 

k Boltzmann's constant 

l mean free path 

m particle mass 

n particle number density 

p, p° local and equilibrium pressures, respectively 

o. 

P traceless pressure tensor 

Np r Prandtl number 

q, q v conductive and convective heat flux, respectively 

q r radial component of q 

Q total heat flux 

Q r radial component of Q 

r radial vector ix + jy; magnitude r is radial distance from pipe axis 

R rate of heat transfer per unit area to pipe 

s entropy density 

s(°) equilibrium entropy density 

s c entropy source strength (collisions! production rate of entropy density) 
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t 


time 


T,T° 

u 

if 


v 


V 


W 


x,y,z 


x i’ x j’ x k 


X 

a 


y 

5: 

e 

■n 

9 

X 

P 


local and equilibrium temperatures, respectively 

reduced particle velocity relative to mean gas flow (see eq. (17)) 

unit tensor 

local flow velocity 

mean particle speed at equilibrium 

magnitude of flow velocity at axis of pipe (r = 0) 

scalar defined by equation (50) 

Cartesian coordinates; also used to denote vector and tensor components 

same as x, y, and z, respectively 

body force per unit mass 

flow parameter defined by equation (35) 

reduced flow velocity defined by equation (33) 

reduced particle velocity relative to mean gas flow; see equation (27) 

Kronecker delta function 

azimuthal angle in particle collisions 

viscosity 

angle between z-axis and particle velocity before collision 
thermal conductivity 
mass density 
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r collision time 

<p first-order (or second -approximation) perturbation function 

<f>,$* dissipation functions 

X scattering or deflection angle 

Special notations: 

collisional time rate of change of f 

0 

substantial derivative - that is, — + v • V 

’ at 

velocity average using f 
velocity average using f(°) 

Primed symbols in collision integrals refer to quantities after a collision; 
unprimed symbols to quantities before a collision. Magnitudes of vector quantities are 
denoted by the same symbol without the arrow. 

THE COLLISION MODEL 

The purpose of this section is to develop a modification of the Krook collision 
model which can serve as the basis of the approximate kinetic equation employed in the 
present research. In particular, the real gas is replaced with a hypothetical system 
having the same gross collision properties but which simplifies the kinetic theory to per- 
mit exact solvability for the velocity distribution function. The reliability of this pro- 
cedure is tested in subsequent sections by calculations of the Prandtl number and com- 
parisons of flux-force relations with Grad’s 13 -moment results. 

A convenient hypothetical gas system for reducing the kinetic equation to a solvable 
differential form is constructed as follows: Introduce a fictitious second gas of infinitely 
massive particles at rest relative to the mean velocity of the real gas at every point; 
assume collisions between real -gas particles to have negligible macroscopic effects com- 
pared with real-gas - fictitious -gas interactions. The collision integrals in the corre- 
sponding Boltzmann kinetic equation are especially simple and resemble in many respects 
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the Lorentz collision integrals employed for certain types of plasma. (See refs. 8 to 11.) 
They comprise the complete collision model used in the present research. 

The utility of the hypothetical gas depends upon its macroscopic correspondence to 
the real gas alone. Although the detailed collision processes and the values of individual 
transport coefficients are different in the two systems, the same general transport phe- 
nomena must occur in both. Undesirable phenomena such as binary and thermal diffusion 
of the real gas relative to the fictitious gas are prohibited by the restriction to the same 
local flow velocity. Hence, the substitution of the values of real-gas transport coeffi- 
cients into the final exact macroscopic equations derived for the hypothetical system 
should yield reliable expressions for the real gas alone; in particular, valuable informa- 
tion should be obtained on the nature and importance of previously neglected contributions 
to flux-force relations from higher Grad moment and Chapman-Enskog approximations. 
The one basic assumption in this statement is that the collision model affects only the 
values of the effective collision frequencies and not the relative importance of transport 
terms, a critical test being the accurate prediction of the Prandtl number. Such numbers 
(that is, ratios of collision frequencies) are generally constant over a large range of 
mass differences between constituent species. 

A simple but informative preliminary illustration of the mathematical character of 
the collision model follows from the application of temperature and pressure gradients 
to the hypothetical gas in such a way as to cause no average flow with respect to the 
laboratory frame of reference. With the concept of order as defined in the "Introduction" 
(that is, with the first-order perturbation function cp being proportional in this case to 
the applied pressure gradient), the corresponding first-order kinetic equation for real- 
gas particles can be written in the following manner for steady-state conditions and no 
body forces: 


Z • Vf = Z • Vf(°) + f(°) c • V<p = 


( 1 ) 


where <T is the particle velocity, is the collisional time rate of change of the 

velocity distribution function f, and 


f = f(°)(l + <p) = n( — 2L_\ 3/2 
\27rkT/ 



<P) 


( 2 ) 


If the temperature and pressure are assumed to vary only in the z-direction and to 
have no second or higher spatial derivatives, the gradient of (p can be deleted from 
equation (1) to yield 
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8f(«) C 2 f<°) 


dz 


( me 2 5W , J_ 9£| = fdf) 
\2kT " 2/9z nk8z] Vat/ C 


( 3 ) 


with the aid of the ideal gas law p = nkT. Finally, the collision term becomes 

= - y [ff (fictitious) - fT (fictitious)] |c - c (fictitious)] b db de dc (fictitious) 

= - J* (f - f T )f(fictitious)cb db de dc(fictitious) 

= -n (f - f’)cb db de = -nef ^ (<f> - 0')t> db de (4) 


according to the properties previously outlined for the fictitious gas and its interactions 
with particles of the real species. Both the real and the fictitious species are assumed 
to have the same number density n in order to facilitate the ultimate substitution of 
the real system for the hypothetical one. 

A convenient trial solution of the resulting kinetic equation 


is 



(<p - cp')b db de 


( 5 ) 


<P = g(c) c z 


( 6 ) 


where g is an unknown function of the magnitude of c and thus is not changed by colli- 
sions with the fixed scattering centers of the infinitely massive fictitious gas. Accord- 
ingly, equation (5) becomes 



As shown in references 8 to 11, the Lorentz-like collision integral in equation (7) 
satisfies 


I (c z - ci)b db d 6 = 


( 8 ) 


8 


if the interaction potential corresponds to Maxwellian particles — that is, if the force 
between particles varies as the inverse fifth power of the separation. This requirement 
is a necessary condition (ref. 8) for the Krook model to describe precisely the collision 
dynamics of the hypothetical system, but it does not diminish the scope of the present 
research because again only the values of transport coefficients are ultimately affected. 
The collision time t in equation (8) is appropriate for the relative diffusion of the two 
species if such diffusion were permitted to occur. 

The substitution of equation (8) into equation (7) yields 


Thus, 



from equation (6). In addition, 


( 9 ) 


( 10 ) 



(U) 


is obtained from equations (4), (8), and (10) and identifies the present method with the 
Krook collision model in this problem. 

Krook' s value of the thermal conductivity is computed from equations (2) and (10) 
and the use of Fourier's law of thermal conduction in the following manner: 

_ 5kpr/9T _1_ 9j)\ 9T _ 5pt 9p q, 

2m \9z nk 9z/ 9z 2nm 9z ' * 


where 


\ = 


5kpr 

2m 


(13) 


Equation (13) differs from the standard expression (ref. 1) by a factor of 2“^/^. Although 
an exact conformity can be obtained with a redefinition of r, this redefinition is not 
essential for present purposes. 
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The complete set of collision integrals comprising the definition of the collision 
model in problems of steady-state viscous flow driven by an applied pressure gradient 
requires a consideration of the following first-order kinetic equation: 


c • Vf(°) + f<°)c • V0 = -njc-v|f(°) J (0 - <p')b db de (14) 

Precise details clearly depend on the exact nature of $ and are deferred in part to the 
next section; however, a significant amount of information can be obtained from the 
Grad 13 -moment velocity distribution function (ref. 9) 


f = f( 0 \l + 


0) = nf-S-) 372 e 
\27ikT/ 


•u 2 1 + |-^ u ^ - l^u^if) • u + i- P : u u 


where P is the traceless pressure tensor defined by 


P = 2p (< u u> 



and 



(15) 


(16) 


(17) 


is a dimensionless particle velocity relative to the mean gas motion. 
In particular, the integral 


J (uu - u' u') b db de 
along with equation (8) written in the form 

JCu - u')b db de = -i-(J!L) 1/2 Ti (18) 

are important in viscous flow. This set is the complete model only if the 13-moment 
approximation is used. 

The evaluation of the new collision integral is very simple in the case of the hypo- 
thetical real-gas — fictitious -gas system. Since U 2 is related to u z by the expression 
(ref. 8) 

u' z = U z cos X - U sin x sin 0 sin e (19) 
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there results 


J( u z - u z 2 ) b db de = + ( u z “ u2 ) sin 2 ej (l - cos2 x ) 

+ 2uu z sin x cos X sin 6 si 11 b db de 
= tt^3u z - u^Jj (l - cos^ x )b db 

= 7tA2 ^3u z - u^* (1 - cos x)b db 

= ||( 3U i ' U2 )j K " U U b tlb * 

' ; fe) 1 / 2 K-“ 2 ) 


2nru\ 


( 20 ) 


with the aid of equation (18) and the definition 


Ai = 


^ (l - cos^x)t> db 
- cos x)t> db 


( 21 ) 


The remaining tensor elements are derived in similar fashion to yield 
jfau - uT!')b db dfi = ^;(i) 1/2 (l35 - u2u) 


( 22 ) 


A preliminary evaluation of the Prandtl number is made by comparing equation (18) 
with the xz-element of equation (22) for a viscous flow in the z-direction. The principal 
difference between the two collision integrals lies in the values of the coefficients; hence, 
the effective collision time for the rate of change of u x u z appears to be 2 /( 3 A 2 ) times 
that of u. Also, since these integrals are the collision integrals associated with the coef- 
ficients of thermal conductivity and viscosity, the same ratio 2 ^ 3 A 2 ^ is expected to be 
the value of the Prandtl number for the hypothetical real-gas — fictitious -gas combination. 

The exact solution of equation (14) merely confirms this prediction and, together with the 
o 

computation of P xz , is deferred to a subsequent section. 
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Meanwhile, a predicted Prandtl number of 2/(3A2j is observed to be very close 
to the Chapman-Enskog value of 2/3 because the calculated value of A 2 is 1.034 for 
Maxwellian particles (ref. 9). Since the Chapman-Enskog result refers to the real -gas 
system and thus involves the rigorous Boltzmann integrals for collisions between real- 
gas particles, an important justification is obtained for the basic assumption of the pres- 
ent research - namely, that the collision model does not affect the relative importance 
of transport terms. The Krook model, on the other hand, yields 

J ( u x u z " W) b db de = -i-( ^L) 1/2 UxU z (23) 

from equations (11) and (15) and therefore corresponds to a Prandtl number of unity. 

The only remaining collision integral associated with the exact solution of equa- 
tion (14) is shown in the next section to be 

I( u z - "i 3 ) b M * = S^fe) 1/2 (( 6A 3 - 3 )" 3 - 3 ( A 3 • 1 ) A Z ] (24) 

the evalutation of which is analogous to the treatment used in equations (19) and (20). 

Equations (18), (22), and (24), together with the statement that u equals u’, thus 
contain all the information necessary for a complete description of the collision model 
employed in the present research. The modifications to the Krook model are systemati- 
cally found from exact evaluations of Lorentz-like collision integrals and correspond to 
the introduction of different effective collision times for the transport of different macro- 
scopic gas properties. A convenient summary of these modifications is provided in the 
13-moment approximation by the following expressions for the Krook collision model and 
the present collision model, respectively: 

(l) c = -7( f - f<0> ) <25) 

and 

(i) c = “?( f ' f (0> ) ~ 2pr f(0) - u2 U) - 2uu] (26) 

Equation (26) is derived from equations (15), (18), and (22). 
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THE EXACT FIRST-ORDER SOLUTION 


A convenient form of the first-order kinetic equation for viscous pipe flow is 
obtained from equation (14) by introducing the equilibrium temperature T° into the 
dimensionless particle velocity 



( 27 ) 


and requiring the perturbation function $ to satisfy 


This procedure permits the spatial dependence of T for a given velocity profile to be 
determined from the solution of the kinetic equation instead of being regarded as input 
data. 


The number density n can be treated in similar fashion, but such treatment is of 
no advantage in the present problem because n is independent of the axial coordinate z. 
(See sketch (a).) This result is required by the equation of continuity 

+ V • (riv) = 0 (29) 

if the system is in a steady state and if 


v = kv(r) 


(30) 


where r is the radial distance from the pipe axis. Moreover, if the induced radial 
variations of n and T are relegated to second order (that is, proportional to the 
square of the flow velocity as partly implied by the energy equation and Fourier’s law), 
the number density is constant in equation (28). 


<my//////////// ///////////////////////^ 


2a 


77777777777777777777 / 77777 / 7777777777777 / 

Sketch (a).- Diagram of viscous pipe flow. The dashed curve 
represents the velocity profile. 
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These restrictions to a steady-state condition, equation (30), and constant n 
through first order, together with the assumptions of no body forces (for example, 
gravitation) and a parabolic velocity profile 


v = 



( 31 ) 


where v 0 is the value of v on the axis and a is the pipe radius, constitute a com- 
plete set of input data for the solution of equation (14) by equation (28). Whether such 
conditions are possible in principle is left for the kinetic theory to answer (by permitting 
or rejecting a solution), as are the predictions of T - T°, transport coefficients, and 
flux-force relations. Velocity and temperature jumps at the walls are ignored in this 
treatment, but can be inserted directly into the final macroscopic expressions. 

The substitution of equations (27) and (28) into equation (14) yields 


1 [ 8f(°) 

f (o) v x 9 x +y y 


8f(°) \ 
9y I 


+ ~y * V(f> = 


i ( y art 

\ r x 8x y Y dyj 


y +y . V<P 


= -n y 


j* (<P - 


db de 


(32) 


if terms containing v 2 are neglected and if ^ is the magnitude of the dimensionless 
flow velocity defined by 



Equation (32) can be written in the more convenient form 


(33) 


2a^xy x + yyyj y z + y • V0 = -n y J (<p - (p')b db de 


(34) 


by using the following combination of equations (31) and (33): 


2£ = -^/_m_\ 1/2 x = c « x 

9x a 2 \2kT 0 / 


(35) 


Since the collision integrals of equations (18), (22), and (24) can be adapted to the 
present notation merely by substituting T° for T and y for u, a simple inspection 
of these expressions and their role in balancing the left-hand side of equation (34) yields 
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i>= (gi + r \)y z + ( x >x + yy y ) % + + + ( x >x + y^e 


+ g 7 r^ + zg 8 r z 



(36) 


as a trial solution. The unknown spatially independent factors g^ are functions only of 
the magnitude of y and thus are not affected by collisions according to the principles of 
energy and momentum conservation in the present model. 

The substitution of equation (36) into equation (34) then gives 


(vfV 


g 2 )( x r x + y y y )y z + (r 2 - rl)s z r z + 2z(xy x + yy y )g 5 


+ r %y z 


+ z (y 2 - r 3 ) g 6 + (xy x + yr y )g 6 r z + g 7 r 3 + 2zg 8 y2 + 3 Z 2g 9 y z J 

= -(gj + r2g 2 )y z - ^^(xy x + y y y )y z - - z ) y l ' 3 ( A 3 - l )v 2y z] 

- z ( xy x + yy y )H - — - ^ 2 ) - zJ Wz 


(37) 


so that 


* ' ( 2 “ (V ' H ‘ + 10ri ] 82 ) 

♦ ^ (2^) V2 g2 . ^(^) 1/2 (» + 2g 2 )(xy x ♦ y , y)rz 


+ (r 2 - 2z 2 )g 2 r z + 2zg 2 (xr x + yr y ) - fr(~j 1/2 (3r 2 - 2z 2 )g 2 


(38) 


from the equating of coefficients of like powers of the independent phase-space variables 
in equation (37) and the use again of equation (36). 

Equations (28) and (38) represent the exact first-order (that is, linear v) solution 
of the viscous-flow problem in question - subject, of course, to the limitations of the 
approximate collision model. This solution is more complete than that of Chapman and 
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Enskog at the same level and for the same model because spatial derivatives of <p are 
retained on the left-hand side of the first-order kinetic equation. 

Conditions on the unknown function g 2 are obtained by taking various velocity 
moments of the distribution function. For example, the relations 



e 


-y2 9 

7 y s 2 


dy = 0 


(39) 


and 



e~y 2 y 4 g 2 dy = 


0 


(40) 


follow from the requirement that the entire contribution to n should come from the 
Maxwellian contribution to equation (28). These integrals further demand zero values 
f° r (Y x )> (^y)’ 311(1 8T / 8r in this problem. 

The consistency of the problem constraints is thus assured if g 2 also satisfies 
the requirement 

<r*> - S jV <* - iff 7 * f e ~ y2y6 ( 2a - 

45it ' A2nm 

= ■ 32 P° t2 ■■ (*! ! - /2 <x _ J°° e ^ 2 y 6 g 2 dyj = 0 (41) 

Agnm V ® 0 / 

that is, if 

e-y 2 y6 g2 dy = ^2 (42) 


Even though equations (39), (40), and (42) do not determine g 2 uniquely, they are 
sufficient to describe completely the usual macroscopic properties of the gas. One illus- 
tration is the local temperature defined by 


T = <y 2 > = T° 


1 + 


32tz 
97T 1 / 2 A 2 


(tt) ' dy 


= T°^l + 


_z_ 9T\ 

ipO Bz J 
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(43) 



so that 


9T _ 32T°t 

dz 




9^/2 H 
4T°Ta/2kT°\ 1//2 8tC>tv o 


3A« 


^2kT° y 


3A 2 a 2 


(44) 


from equations (35), (42), and (43). 

Equations (43) and (44) can be employed in the following manner to express equa- 
tions (28) and (38) in a more convenient form: 

* <o) ■ r vu2 ■ ^ + ^ 


n/-SS_\ 3/2 flf) 3/2 

\27rkT°/ \T / 


exp 


Jl_(e-v) 2 -?^iu2 

2kT° V ' 3A 2 a 2 




A 2 a 2 


3A 2 a 2 



- - 3)1 

1 

3A 2 a 2 


(45) 


Accqrdingly, the complete first-order solution is written as 

f = f(°)(l + <p) 


(46) 


where f(°) is defined now by the first equality in equation (45) and the new <f> satisfies 
the expression 


<P = 


8 tv^ 


3A2a^ 


(" 2 - 1) - SaTT^Hv 2 - "IK * (*”* * 


_ 3 \ . 

2 T 

2/ 

5A 3 - 

4 tzu 2 

/2kT\ 

*2 

\^r) 


8prV 


z \ 3A2(5A 3 - 3)nm| 


| 5A 3 - 9)u 2 


+ 10u |s 2 + -^dr) V2 B2 - ^(^5 t) 1/2 ( xu x + y u y)g 2 u z + ( r2 ' 2z2 )g 2 u Z + 2z S 2 ( xu x + 3™y) 




( 47 ) 
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This formulation of the velocity distribution function clearly shows an important 
distinction between the Chapman -Enskog first-order perturbation function and that of 
equation (47). Whereas the Chapman -Enskog method has the primitive property of rele- 
gating the entire spatial dependence to the first five moments n, T, and v, numerous 
terms in equation (47) involve x, y, and z explicitly. Hence, many higher moments 
must be introduced and permitted their own spatial variations before the first-order dis- 
tribution function is completely described by local values of velocity moments. The 
remainder of the present section is concerned with several steps toward such a 
representation. 

The next higher vector moment is the conductive heat flux given by 

i . - *) 2 (* - <^> ■ ki ?) i/2 i “ 2 « * 

= iKrar) 1 ^ 1 Z?+ ^ ‘ 2Z % e “ u2u6 S2 du 


£ 28p2r\ r ga 2 

3A 2 nma 2 _ 


21v, 


4 


2kT^/ 2 

7rm / 


£° e“ u2 u 8 g 2 du 


P v. 


2z? + k[r2-2z 2 + |?EI 


3A2nm i 


- £ I^V/atTV 172 r e -u 2 u 8 g du (48) 

9A 2 nm\7rm / Jq & 


Besides the direct use in this expression of equation (47) to evaluate the original integral, 
equations (42) and (44) are employed in the subsequent simplifications. 

If q and its spatial derivatives from equation (48) and the spatial derivatives 
of v from equation (31) are substituted for certain combinations of terms in equa- 
tion (47), there results 


<P = - 


8V 2 hx "tfl fav 2 (fv is5|\ fi . a2f2kT\V 2 1 
ax + 5p \az + ax / ^ + ay 5p\az ay /j^ [_ v o ' m ' 2 2 


20t i 


27A 2 ^ 

ax 5p\az 

4t 1 

f u 2 . 3 u 2^ 

15A 2 p 

l u dz 


^/_m-\!/ 2 / 2 _%.-5 
5p\2kT/ \ 2 r 


16pr 2 v 0 u z 


15A2(5Ag - 3)nma‘ 


:(if) ( 4 ( 5A 3- 3 )( u2 -f) 

(Equation continued on next page) 
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as a partially complete moment representation of the perturbation function. 

Another step toward the complete moment representation of the first-order distri- 
bution function is suggested by the last term in equation (49). Since the components 
of u do not appear in this term, it is convenient to introduce the scalar W according 
to the definition 


W = <u 4 > - <u 4 ) 0 = i j* u 4 (f - f<°))dc = i J u 4 f(°4 dc (50) 

This definition makes W the next higher scalar moment beyond (u°) corresponding 
to the number density n and (u^) corresponding to the temperature T; hence, the 
mathematical sequence is well defined even though the moments higher than are 

not familiar physical entities. 

The evaluation of W by the use of equation (49) in equation (50) and the aid of 
equations (44) and (48) yields 
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This expression and 



Despite the fact that the explicit spatial dependence has been removed in equa- 
tion (53), that expression is still somewhat complicated because the function g 2 is 
unknown. The conditions given in equations (39), (40), and (42) do not determine g 2 
uniquely; thus, the present problem is not completely defined by the restrictions to a 
steady state, constant number density, and parabolic velocity profile. Fortunately, how- 
ever, the choice of g 2 has almost no practical significance because of the following 
summary of the effects of different g 2 selections on the usual 13 moments: none at all 

on n, T, v, and P if equations (39), (40), and (42) hold; almost none on q as is shown 
in a subsequent section. 
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is 


The simplest form of g 2 which is consistent with equations (39), (40), and (42) 




V2, 


Accordingly, since 


fe-uV 




du = 


63v, 


8a 2 


o /77-m \ 

2 \2kT/ 
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equation (53) becomes 
8tu„ 
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. 5u 2 + 25\/fiv + + + [av + 2_/^jr + ^ 
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15nm 


lSA^p 


+ -r- — llu^ - 45u 2 + - 6^u 4 - 5u 2 + 


9z 


4t 2 v 0 u z 


( 54 ) 


(55) 


(56) 


e / — |(5 A 3 - 3) (llu 4 - 43u 2 + 

15A 2 (5A 3 - 3)a 2 '' m / \> 3 / V 4/ 

- 2 [( 5 A 3 - »)u 2 + 10u|](u« - 5u 2 + 

THE TRACELESS PRESSURE TENSOR 


The traceless pressure tensor defined in equation (16) is computed by equations (44), 
(48), and (53) to satisfy 




(Equation continued on next page) 
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if the viscosity 77 is given by 


(59) 

and if proper account is taken of the characteristics of v in the present problem. 

This value of the viscosity combines with the thermal conductivity of equation (13) 
to yield the correct Prandtl number 


N - C P ?? - 5kr? 2 2 

Pr X 2mX 3A 2 ~ 3 


(60) 


thereby at least one principal improvement of the present collision model over that of 
Krook is confirmed. Equations (57), (58), and (60) thus provide some support for the 
hypothesis that precise flux-force relations and ratios of transport coefficients for real 
gases can be obtained from the kinetic theory of the real-gas— fictitious-gas combination. 

Although the present formulation of the traceless pressure tensor is the same as 
the 13 -moment approximation, the evaluations of the tensor elements are different because 
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of different results for the spatial derivatives of q. This aspect is examined in some 
detail in the next section, where a more convenient representation of q is deduced. 

CONDUCTIVE HEAT FLUX 


An expression for the conductive heat flux q has already been derived in equa- 
tion (48), the combination of which with equations (13), (44), and (52) gives 


5 - g 4 p 2 r 2 v 0 2 p 2 T 


A2nma^ 


3nm 


VW = -k ^ — - ^3- VW 


5 8z 


15 


= -kx2T + £ »8T.4^r w 
8z 5 dz 15 


(61) 


Since the divergence of P satisfies 


o 

V • P = k| 


(dp. 


xz 


9P 


TL 


o 

8P r 


zz 


= -knk 9T = -nkVT 
8z 


^ 8x ' 9y 9z 
from equations (31), (44), (48), (58), and (59), equation (61) becomes 


(62) 


q = -AVT - V • P - ^7? W 
5nk 15 


(63) 


A convenient comparison with the Grad 13-moment result is obtained by taking the 
mc 2 cT-moment of equation (1) and using either equation (15) or equation (56) in the present 
collision model. The result is 

V • (nm(c^cc)) = V • j ^ n ^ T ((u 2 uu) + /3 (u 2 u) + {u 2 u)js + 2 /3 • (UiTu ) jj 

- 4 # V ' (t 2 <u%U>) . £ 8#£<*j> 0 <u2uu> 

= T VT + HnT V ' ( u2 ““> = m j" d? 

= — ^2 — f e" u2 u 3 u($ - 0 f )b db de du= -^q (64) 

7T 3 / 2 m J 

Second-order terms are neglected in the simplifications of V • (nm(c 2 c"<f )), whereas 
equations (4) and (18) are employed in the evaluation of the collision integral. 


23 



Equation (64) yields 


q = -a&|L VT + V • <u 2 uu) q + V • (<u 2 uu> - <u 2 uu> 0 )J 
= -2XVT - ^ v • (<u 2 uu> - <u 2 uu) q ) 


(65) 


with the aid of equation (13). Unlike the collision integral in equation (64), the value of 
which is independent of whether equation (15) or equation (56) is used, the evaluations of 
(u 2 uu) - (u 2 uu) n are 


{u 2 uu^ - (u 2 uu) = — P 
° 4p 


( 66 ) 


and 


(u 2 uu) - (u 2 uu). = — — fl7 P + lofiiP 
° 36p|_ \ 


Xa O aaO 


P xx + U p yy + kkr U, 




(67) 


respectively, for the Grad 13-moment and the present velocity distribution functions. 
Accordingly, the two determinations of the conductive heat flux are given by 

q(13 moment) = -AVT - — V • p" 

5p 

and 


( 68 ) 


q(exact first order) = -AVT - V • P - VW 


(69) 


from equations (44), (48), (58), (62), and (65) to (67). Equation (68) is also identical with 
the Grad 13-moment result quoted in reference 1 (p. 495; note, however, the typographical 
error). 

The fact that equation (69) reduces to equation (68) in the 13-moment limit of W = 0 
completes the theoretical support for the hypothesis that precise flux -force relations for 
real gases can be obtained from the kinetic theory of the real-gas — fictitious -gas combi- 
nation. Of the 13 physically important moments n, P, q, and 




( 70 ) 
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from equations (31) and (44), the formal results of the present research differ from those 
of the Grad 13 -moment approximation only in the expression for q; this one difference 
stems not from the collision model, but from the inaccurate 13-moment determination of 
the quantity (u 2 uu) in equation (65). 

Consequently, the six moments corresponding to (u^uCT) and three more to cover 
the u • VW term in equation (56) are minimal additions to the 13 -moment distribution 
function in order to provide an adequate first-order description of the present problem. 
The flux -force relations of equations (58), (69), and (70) are thus equivalent to a Grad 
22 -moment approximation or better. 

This discussion would be trivial, of course, if the last term on the right side of 
equation (69) were negligible compared with the first two. Actually, the reverse is true 
because of the following combinations (correct through first order) of equations (13), (52), 
(55), (59), (62), (69), and (70): 

VW = 15pV . ^ |2zr + k(r 2 - 2z 2 )l - V • 1? (71) 

4XTa 2 L J 4 P 


and 


q = -XVT + V 
5p 


f P v ( 
P 


pj2zr + k(r 2 - 2z 2 jj 


= - 2z, - 5 gF) 

Not only does the VW contribution to q add a term of ordinary conduction magnitude 
(that is, XVT), it also adds much larger terms of convection magnitude which are inde- 
pendent of the choice of g 2 in the perturbation function. An interesting curiosity of 
equation (72) with regard to the second law of thermodynamics is given in the appendix. 

Since the convective heat flux is deduced from equation (70) to be 






5pv 0 (a 2 - r 2 ) . 

« k 

2a 2 


(73) 


the total heat flux Q can be written as 


Q = q + q v * 



(74) 
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Consequently, the ratio of the maximum r-component to the maximum z-component of 
Q is 

Qrfr = ») 4az „„ 

Q z ( r = 0) 5a^ + 4z^ 

Except at z = 0, which corresponds to T = T° from equation (43), the picture of 
heat transfer obtained from equations (72) to (75) is entirely different from the 13 -moment 
result of Q ~ q y . This difference is especially noticeable in the first-order expressions 

Q r (13 moment) = 0 (76) 


and 

Q r (exact solution) = - 2p J° rZ = p2 ^ T j T T ^ C 77 ) 

from equations (43), (70), and (74). Whereas the 13-moment approximation predicts no 
heat transfer to the walls of convection magnitude (that is, first order or proportional 
to v), the present solution predicts an outflow of energy in the heated regions of the gas 
and an inflow of energy in the cooled regions - even though T (first order) does not 
vary with r. 

It is important also to realize that the prediction of a first-order radial heat trans- 
fer has nothing to do with the usual macroscopic energy equation, which is concerned with 
the divergence of q and yields a second-order or v^ -contribution to q r . Equation (72) 
states that the spatial variations of q_ and q_ are such as to make the divergence 
of q vanish in first order; hence, the phenomenon discovered here is indeed new in the 
sense that it cannot be predicted for this specific problem by either the Navier-Stokes 
theory or the Grad 13-moment approximation. The latter method modifies only the neg- 
ligible contributions to Q. 

As mentioned previously, equation (72) also has important effects on the values of 
the traceless pressure tensor elements in equation (58). Whereas the 13-moment approx- 

imation to P satisfies the Navier-Stokes expression 

-,[t£ + £i)g + (j£.£f)|| 

= ^Sp|(ik + ki)x + ( jk + kj)y] (78) 
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in the present problem because of the vanishing of the spatial derivatives of q 
(13 moment) to first order, the exact solution gives 



Accordingly, major changes are made in Newton's law of viscosity as well as in 

Fourier's law of thermal conduction; in particular, numerical changes occur in the off- 

o 

diagonal elements of P and nonzero diagonal elements are introduced. However, no 
change appears in the important relation 

-£L - 4tjv 0 , , 

V • P = k — ^ = -nkVT (80) 

a 2 


which is fixed by the macroscopic equation of motion 


p ^ + v • vjv - pX + Vp + V • P = Vp + V 


P = 


( 81 ) 


THE ENERGY EQUATION 


The final computations of the present research concern the energy equation, which 
is obtained as follows by taking the mc 2 -moment of the kinetic equation and using the 
equations of continuity and motion: 


A(nm(c*0) - 2nmX * v + V • (nm(c2c)) 

= ^jnrn^( c - v)^^j + ~(nmv^) - 2nmX - v + V - jnm^(c - v)^(c - v) + (c - v)^v + 2v • (c - v) (c - v) + v^v^J 


= 3 2E + 2pv • - v • Vvj - v2v • (pv) - 2 pv - X + 2V • q + 3pV - v + 3v • Vp + 2V • jv • (p + pu) 

= 2pc p - 2 gE+ 2V - q - 2^ • (v • f) - V • (v • f)| 


+ v 2 V • (pv) + 2 pv • (v • Vv ) 


( 82 ) 
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Thus 


pc Pm = m ' v '^ +$+$ * (83) 


where $ and $* are dissipation functions satisfying 

$ + $* = ? . (v • p) - V • (v • "P ) (84) 



8vj 8vj 
ij 9Xj + 9Xj 


(87) 


and 


* 9( li 9 <lj 

e ii = 7 — + — “ 

13 8x j 9x i 


( 88 ) 


Third-order terms proportional to v 2 and containing spatial derivatives of t? are 
neglected in equations (85) and (86). 

Except for the dissipation function $*, which derives from the heat-flux coupling 
terms in equation (58), equation (83) is the Navier -Stokes energy equation. The values 
of $ and $* for the two distribution functions are obtained from equations (68), (70), 
(72), and (85) to (88) to be 


$(13 moment) = $(exact solution) = r] 



4?7v 0 2 r 2 

a 4 


(89) 


$*(13 moment) = 0 


(90) 
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and 


$*(exact solution) = 


16?7V 0 2 r 2 

5a 4 


(91) 


through second order. 

The resulting energy equations for the present problem are for the 13 -moment 
solution: 

V • q = -| v • Vp + $ + 4>* = — (a 2 - |r 2 ) (92) 

and for the exact solution: 

V • q = -|- v • Vp + <f> + <3>* = - - ^a 2 + ir 2 j (93) 

Equation (70) is used in expressing Vp in terms of kv Q . 

As in the evaluation of the traceless pressure tensor elements, and even though the 
formal expression of the energy equation is the same as the 13-moment approximation, 
the new result for q given in equation (72) has pronounced effects on the values of <$* 
and the divergence of q. Equation (72) and the solutions of equations (92) and (93) imply 
not only the substantial changes discussed previously with regard to the Navier-Stokes 
picture of first-order heat flux, but also significant alterations in the v 2 -dependent con- 
tribution to q r . 

The rate of heat transfer per unit area to the walls of the pipe is defined by 


R = - — / 27raQ r (lst order; r = a)Az + \ V • Q(2d order)dx dy dz 

27ra Az J 

- 27 jQ z (2d order; z + Az) - Q z (2d order; z)Jr drj 

i r a — 

~ Q r (lst order; r = a) + — \ rV • q dr 

a Jq 


(94) 


where Az is a small increment of axial length. Two assumptions are made in the 
simplifications of equation (94): (1) the convective contribution to V • Q (2d order) is 

negligible or zero, and (2) changes in Q z (2d order) along the z-axis are third order. 
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Both are reasonable in view of the fact that the flow velocity has no radial component and 
the likelihood that the nth derivative with respect to z corresponds to a term of 
nth order. 

Equations (43), (70), (74), and (92) to (94) thus yield 


R(13 moment) = 


5t7v 0 2 

2a 


and 


R(exact solution) = - 


2pv 0 z 

a 


+ 


33t ? v 0 2 _ ap 2 (T - T°) t 3377V 0 2 
10a 2r)T + 10a 


(95) 


(96) 


COMPARISONS WITH CHAPMAN -ENSKOG THEORY 

Correlations of the preceding results with the kinetic theory of Chapman and Enskog 
(ref. 2) are somewhat ambiguous because the use of r (or the mean free path) as a per- 
turbation parameter is not valid in the present problem. The reasons are obvious both 
physically and from the exact first-order solution in equation (47). In particular, the 
approach of r (and thus the viscous forces) to zero is inconsistent with the steady-state 
condition unless the applied temperature gradient also vanishes, for otherwise an unbal- 
anced driving force would exist and cause the gas to be constantly accelerated through the 
pipe. The same conclusion is reached by analyzing equation (47), the last three terms of 
which are proportional to r“l, t - 1, and r _ 2 because of the r-1 dependence of g£ 
from equations (54), (59), and (70). Hence, the fundamental perturbation principle of 
Chapman and Enskog does not apply to this case. 

Comparisons can be obtained, however, without utilizing the explicit and implicit 
appearances of r in equation (47); for example, one can simply ignore the perturbation 
principle and look directly at the Chapman-Enskog perturbation equations for a basis of 
correlation. A particularly convenient and systematic scheme is established as follows: 
first assume that no term in equation (47) occurs below the Chapman-Enskog second 
approximation because of the factor v 0 or g2> next assume that V<fi fails to appear 
in the second-approximation kinetic equation because <fi (second approximation) does 
not contain x, y, or z explicitly; then assume that second spatial derivatives of <fi 
are absent from the third-approximation kinetic equation because <fi (third approximation) 
contains no explicit x 2 , y 2 , z 2 , xy, xz, or yz factors and so forth. Each power of 
x, y, or z appearing explicitly in equation (47) thus represents a raising of the 
Chapman-Enskog level of approximation by unity in order to include the corresponding 
term. As explained previously, the occurrence of such higher approximations in the 
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present first-order solution derives from the use of a more complete first-order (or 
second- approximation) kinetic equation which includes all terms linearly proportional 
to the viscous flow velocity. 


The assignment of a level of approximation to a particular form of a flux -force 
relation must also involve the selection of the pertinent terms in equation (47). For 
example, the VW - contribution to equation (69) derives from the (r^ _ 2z^)g2^ z and 
2 zg 2 ^xu x + yuy) terms in equation (47); accordingly, that expression for the conductive 
heat flux is assigned to the fourth Chapman-Enskog approximation. A complete summary 
is given in table I. 


Two exceptions to the correlation scheme are noted in table I, each of which results 
from inconsistent Chapman-Enskog assumptions. Equation (47) indicates that the first 
appearance of a nonzero traceless pressure tensor should occur in the third approxima- 
tion, whereas Chapman and Enskog credit Newton’s law of viscosity to the second approx- 
imation; in addition, equation (47) clearly shows that the coupling terms corresponding to 
the Grad 13-moment and Burnett formulations cannot be a priori separated from Newton's 
contributions. Although spatial derivatives of the viscous flow velocity result from spatial 
derivatives of Chapman and Enskog's Maxwellian distribution function in the second- 
approximation kinetic equation, they should be deleted after the fact as being inconsistent 
with that equation. The reason is obvious: Since the resulting <£ (second approximation) 
must contain a term proportional to 3v/9x, which, in turn, is proportional to x itself in 
the case of a parabolic velocity profile, the neglect of 30/3x in the second -approximation 
kinetic equation is inconsistent with the assignment of Newton's law of viscosity to that 


TABLE I.- SUMMARY OF MACROSCOPIC RELATIONS 


Number of moments 
in Grad approximation 

Level of Chapman- 
Enskog approximation 

Pressure tegsor 
element, 

Heat flux, q 

Energy equation 

5 

First 

0 

o 





.... 

Second 

/9vi 9vj 2 * „ -A 

-XVT 


13 

Third 

'^9Xj 8x t 3 l J 7 

5p\3xj 8Xj 3 U J 

-XVT - V * ~P 

5p 

pc Pi?=m- v - g+ * + *‘ 

*§22 

Fourth 

2« v-v) 

y3xj dx i 3 y / 

-XVT - V • *P 

5p 

pc p^ = m ■ v ^ + ® +$ * 



2 

5pyaxj dxj 3 J 

-~W 

15 



*The value of the VW term in the equation for q depends upon the number of moments employed. At least 
22 moments are required for this term to appear. 
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level. Either the first spatial derivatives of <f> should be included in the second- 

o. 

approximation kinetic equation, in which case the complete 13 -moment formula for P 

would be second approximation, or no contribution to P should occur below the third 
approximation. 

A similar argument applies to Fourier’s law of thermal conduction in Chapman and 
Enskog's second approximation. Equation (47) again implies that the 13-moment coupling 
term cannot be a priori separated from Fourier’s law, both contributions being assigned 
to the second approximation by the present correlation scheme. Chapman and Enskog’s 
neglect of the divergence of the traceless pressure tensor in the second approximation is 
therefore inconsistent in the present problem. 

Not only are the arguments similar, but the two Chapman-Enskog inconsistencies 

„ g. 

with respect to P and q are, in fact, related to each other. If the gradient of (p 
(and thus the gradient of the spatial derivatives of the viscous flow velocity) had been 

included in the second-approximation kinetic equation in order to be consistent with the 

q 

retention of 3v/9x, etc., the divergence of P would have occurred automatically in the 
second approximation to q. 

Completely consistent treatments, whether based on the present correlation scheme 
or something different, always predict more complicated expressions than those of 
Newton and Fourier for the present problem. This conclusion further explains the 
absence in table 1 of a direct correspondence between the Chapman-Enskog second approx- 
imation and one of the Grad N-moment variety, the latter method being free of basic incon- 
sistencies. Only the use of insufficient moments seems to affect the accuracy of the gen- 
eral Grad approximation for neutral gases, although other difficulties arise in certain 
plasma applications. (See ref. 11.) 

CONCLUDING REMARKS 

Exact first-order solutions of an approximate kinetic equation have been found for 
the problem of viscous gases flowing through circular pipes with constant temperature 
gradients. The assumption that errors inherent in the collision model affect only the 
values of transport coefficients and not their ratios or the formal flux-force relations is 
supported in two ways: (1) the expressions for all 13 physically important moments 
reduce in the proper limit to those of the Grad 13-moment approximation, and (2) the 
ratio of the calculated viscosity and thermal conductivity agrees with Chapman and 
Enskog's result. 
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A significant new contribution to the conductive heat flux was discovered which has 
the same order of magnitude as heat convection and substantially changes the Navier- 
Stokes description of energy transfer in the problem considered; in addition, the evalua- 
tions of the traceless pressure tensor and the dissipation terms in the macroscopic 
energy equation are affected. The description of these effects requires at least the 
fourth Chapman-Enskog approximation or the use of 22 moments in the Grad distribution 
function. 

Although the problem considered in the present research is possible in principle, 
experimental difficulties may prevent a direct laboratory confirmation of the results. 

For one thing, a two-dimensional control of the temperature variations within the pipe 
and its environment would have to be maintained in order to guarantee a constant axial 
temperature gradient and also a steady radial flow of heat into or out of the gas. 

In any event, the principal value of the research probably lies more in the method 
introduced than in the actual problem solved. The extension of the collision model to 
more complicated integrals is precisely defined, so that the only question is whether the 
resulting differential equations can be solved for given situations. Since the technique 
of solution is fairly straightforward, no great difficulty is anticipated — even for time- 
dependent systems. The most important precaution is that the first-order kinetic equa- 
tion be complete in the first power of some small physical parameter (the viscous flow 
velocity in the present problem), for that is the device by which significant contributions 
from higher Chapman-Enskog or Grad moment approximations are included with minimum 
effort. This saving of effort is generally very substantial. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., April 8, 1970. 



APPENDIX 


SECOND LAW OF THERMODYNAMICS 

The purpose of this appendix is to discuss an interesting curiosity of equation (72) 
with regard to the second law of thermodynamics. Since the dominant contribution to the 
axial component of q is given by 

q. = -E£(r2 - 2 z2) <A1) 

a* 

and since v Q is positive if the temperature decreases with increasing z, conductive 
heat flows from cooler to warmer regions in the vicinity of z = 0. This result directly 
contradicts one of the simplest and most popular interpretations of the second law. 

An explanation of this paradox is obtained from the following more basic and cor- 
rect formulation of the second law: The entropy source strength (that is, the rate of 
entropy production per unit volume) must be nonnegative at every point. The statistical 
expression for the entropy source strength is deduced from the local entropy density 

s = s(°) - k f f log - 7 — dcf (A2) 

J e f (o) 

to be 

= ' k I [f + l0ge ^ + t) d?= ~ k I l0ge ^ + ^(l) ^ ~ " k I ^(ff) c ^ A3) 

If the left-hand side of equation (1) is substituted for ( 8 f/ 8 t) c in equation (A3), the 
result can be written 

s c = -k y (pc • ^Vf^ 0 ) + f (°)v<pj dc 

- - k J ,(0> *{* ' [(“ 2 - |)¥ + V *| + 2 £ U r“ 2 }« 

= • VT + i(l) 2 - f (w) V2 1 *- u2 * ■ ■ * < A4 > 
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APPENDIX - Concluded 


with the aid of equations (16), (48), (70), and (79). Contributions of order v 8 and higher 
are neglected. Since the last term of equation (A4) never appears in the second Chapman- 
Enskog approximation or Navier-Stokes theory, the requirement that s c be nonnegative 
usually restricts q to the direction opposite of VT. 

In the present problem, however, the last term of equation (A4) is the most impor- 
tant part of s c . This conclusion is reached by first introducing the mean free path 
l = vt, where 


V = 2 feT) V2 ( 

\7rm / 

and using equations (44), (59), (70), and (Al) to rewrite equation (A4) in the convenient 
form 


TS, 

nk 


r “ - hS eV * ' ' * 


(A6) 


Since the function g 2 and the constant a are both proportional to r - * by equa- 
tions (42), (44), (59), and (70), the largest contributions to 4> in terms of small values 
of 1 / a are the last three terms in equation (47). Accordingly, 


0 * ( r 2 - 2z 2 )g 2 u z + 2rzg 2 u r - |-g_y /2 (3r 2 - 2z 2 )g 2 (A7) 

and 

" • V<*> « 4rg 2 u r u z + 2zg 2 (u 2 - 3u 2 ) - g 2 [ 2rzu r + ( r2 " 2z2 ) u z] < A8 ) 

so that 

(A9) 

3 IT ' 

at the position r = a, z = 0 where the previously mentioned paradox with the second 
law is most pronounced. 

No disagreement with the basic formulation of the second law of thermodynamics 
occurs, nor can a violation ever exist when the solution of the kinetic equation is exact 
in the present sense. 
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